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We reconsider the 2d model for CuGeOs introduced pre- 
viously (Phys. Rev. Lett. 79, 163 (1997)). Using a computer 
aided perturbation method based on flow equations we ex- 
pand the 1-triplet dispersion up to 10 th order. The expansion 
is provided as a polynom in the model parameters. The latter 
are fixed by fitting the theoretical result to experimental data 
obtained by INS. For a dimerization S ~ 0.08(1) we find an 
excellent agreement with experiment. This value is at least 
2 to 3 times higher than values deduced previously from Id 
chain approaches. For the intrachain frustration «o we find 
a smaller value of 0.25(3). The existence of interchain frus- 
tration conjectured previously is confirmed by the analysis of 
temperature dependent susceptibility. 

75.40.Gb, 75.10.Jm, 75.50.Ee 



I. INTRODUCTION 

The dispersion of the magnetic excitations is an im- 
portant source of information on experimental low- 
dimensional spin systems. Knowledge of the dispersion 
relation u>(k) helps essentially to identify the model ap- 
propriate to describe the compound under study. The 
dispersion relation provides also important insight in the 
nature of the ground state. Very common in low dimen- 
sional systems is the scenario of a singlet S — ground 
state without magnetic long range order (a "spin liquid" ) 
of which the elementary excitations are triplets S = 1. 
These systems are generically gapped. Examples are iso- 
lated or weakly coupled dimerized spin chains and spin 
ladders such as (VO) 2 P 2 7 , § the spin-Peierls phase of 
CuGe0 3 , @,§ and SrCu 2 3 . @ A true 2D example is 
SrCu2(BO)2 which is characterized by frustrated dimers. 

M 

In these gapped S =1/2 systems where the gap is re- 
lated to some "strong" bond (which can be also the rung 
of a 2-leg ladder) the elementary triplet excitations are 
in principle accessible by a perturbative expansion about 
the limit of isolated dimers. This approach, however, be- 
comes tedious for the description of realistic materials 
since the expansion parameter is often not really small. 



Thus one has to compute high orders to achieve quan- 
titative agreement. For this reason various automated 
approaches have been conceived which leave the tedious 
part to computers P O]. 



In the present article, we will apply thepreviously in- 
troduced perturbation by flow equation pd[ to the two- 
dimensional, though anisotropic, system of CuGeOs in 
its dimerized low-temperature phase [Q . Thereby we ex- 
tend the previous analysis (Ref. henceforth cited as 
[I]) considerably. Our starting point remains the same as 
before. The strongest coupling is given by J; the other 
couplings are given relative to J as indicated in Fig. |l| 
(for details see Fig. 1 in [I]). 
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FIG. 1. Dimerization pattern in the dimerized low tem- 
perature phase of CuGeOa. The couplings are denoted rela- 
tive to the strongest coupling J which is set to unity in the 
figure. 



II. METHOD 

The problem to be solved reads 

H = Ho + XH S . (1) 

As in the chain in Ref. jl]J the isolated dimer limit (A = 
at finite fi/\) has an equidistant energy spectrum and 
the perturbation can alter the number of energy quanta 
(here: triplets on the dimers) by 2 at maximum. Hence 
H s can be represented as H s = T_ 2 + T-i + T + T x + T 2 
where Tj stands for the perturbing part changing the 
number of elementary triplets by i. The same formal- 
ism as in Ref. can be used. This formalism maps the 
perturbed Hamiltonian by a continuous unitary trans- 
formation, the so-called flow equation method [Q, to 
an effective Hamiltonian H D f{ which conserves the num- 
ber of energy quanta, i.e. = [H c g,Ho\. The effective 
Hamiltonian has the form 



1 



oo 

h cS = H + Y, xk Yl c (m)T(m) , (2) 

k=l \m\=k,M(m)=0 

where to is a vector of dimension k of which the com- 
ponents are in {±2, ±1,0}; A/ (to) = signifies that the 
sum of the components vanishes which reflects the con- 
servation of the number of energy quanta (triplets). The 
coefficients C(m) are generally valid fractions computed 
in Ref. |ll]] where also further details on the flow equation 
method can be found. 

Since if e ff conserves the triplet number the one-triplet 
sector is particularly easy to solve. Acting on one triplet 
the action of H c g may only consist in shifting the triplet. 
This means that the triplet hops on an effective lattice 
where one site stands for one dimer on the original lattice, 
see Fig. 2 in [I] or Fig. |. 




FIG. 2. Effective lattice on which the triplet hops. We 
calculate the amplitudes hjj for all hopping processes starting 
on (0,0) and ending on one of the depicted dimers (circles). 
For all circles that are accessible by an arbitrary hopping of 
length six or less the amplitudes have been calculated within 
6 th order. (The length of a hop is the minimum number of 
bonds (solid or dashed) required to link start and end point.) 
The amplitudes for light gray, dark gray and black circles have 
been extended in 8 th order, provided that these sites can be 
reached by a hopping oc fj, of length 8. Analogously, the 
amplitudes for dark gray and black circles were extended by 
calculating hopping processes oc jj, 1 within 9th order. Finally, 
the amplitudes for all black circles were extended by processes 
oc fi within 10th order. The arrows indicate axes with respect 
to which reflection symmetry holds. 

The full dispersion oj(k) is obtained by Fourier trans- 
form 

u)(k) = J/^ j hj_ n exp(«(/cij + fc2«)) • (3) 




o o o 
o o 
o 



FIG. 3. Computer generated cluster necessary to compute 
in 8 th order, allowing for arbitrary hopping processes of 
length 8. The light gray (dark gray) circle denotes the start 
(end) dimer (0,0) ((3,-1)). 

The hopping amplitudes hj^ n can be calculated on fi- 
nite clusters of the (in principle infinite) effective lattice: 
From the linked cluster theorem we know that the finite 
order contribution of a short-ranged perturbation does 
not depend on the cluster size for sufficient large clus- 
ters. Carrying out our perturbation within order I im- 
plies that one allows dimer to dimer hopping processes 
of length I p"l| ] . The minimum cluster for a given ampli- 
tude hj^ n in a given order I contains all dimers and links 
that are involved in a hopping of length < I starting at 
dimer (0, 0) and ending at (J, n). The minimum cluster is 
determined by considering all paths from (0,0) to (j,n). 
All dimers and links covered by one of these paths are 
part of the minimum cluster. In Fig. |^, the computer 
generated minimum cluster for calculating /i3,_i in order 
8 is shown. 

Due to the strong anisotropy of the quasi-lD system 
CuGe03 it is reasonable to use higher order terms only 
along the chains. This simplifies the computational task 
considerably since the calculation of a hopping process 
along the chain is much simpler. The cluster to be consid- 
ered can be chosen smaller. The same is true for hopping 
processes close to the chain direction. Here we restrict 
the hopping processes to be at maximum quadratic in 
the interchain hopping /i, which reduces the cluster sizes 
significantly so that the perturbation order can be en- 
larged. 

III. ANALYSIS OF EXPERIMENTAL DATA 

The results for the hj, n are too lengthy to be published 
in written form. We will provide them in electronic form 
on our home pages on appearance of this article. In [I] 
the hj jn in third order in A and fi were presented. A 
few of these are erroneous. They are corrected herewith 
]l3[ . The corrections, however, have no influence on the 
conclusions in [I] (see also discussion below). 

Once all amplitudes hj iTl are calculated the dispersion 
relation is given by Eq. (|^). After rewriting Eq. (||) in 
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terms of fcf, and k c (the reciprocal basis to e& and e c ) we 
add the term At a cos(fc a ) cos(/c c ) with 4t a = 0.22 meV to 
account for the dispersion in a-direction (cf. [I]). To fix 
the parameters J, a, (3, /i and A (cf. Fig. |l|) we use the 
one-magnon dispersion data for CuGeOa experimentally 
determined by inelastic neutron scattering Note 
that the hopping amplitudes are computed as polyno- 
mials over Q in the parameters. 

As noticed in [I] the parameter (3 has almost no in- 
fluence on the shape of the dispersion. Hence we re- 
frain from determining (3 from the dispersion but set 
it beforehand to some reasonable values in the inter- 
val [—0.3,0.3]. This choice is motivated by comparing 
the microscopic direct super-exchange path \x and the 
shifted super-exchange path fi(3 shown schematically in 
Fig. | ( cf. Fig. @). There is only one path per Cu + - 
site for the shifted coupling whereas there are two paths 
for the direct coupling. Thus we expect \fi(3\ w 1/2|//|, 
i.e. |/3 1 « 0.5. There are also results from ab-initio cal- 
culation for the interchain couplings which indicate the 
existence of interchain frustration |lq ] . Further evidence 
is provided below by the analysis of the susceptibility. 
Furthermore, we find that for > 0.4 fits to the disper- 
sion data become worse. 
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FIG. 5. Dispersion ui(k a = 0, fcf, = 0,1/2 + k c ) in 
c*-direction. The arrows indicate the experimental points 
used to fix the parameters. 10 th order fits based on 
(f3,X,a,fi,J) = (0.3, 0.836, 0.225, 0.266, 13.1meV) and (-0.3, 
0.846, 0.209, 0.081, 12.3meV), respectively. 
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FIG. 4. Schematic view of the microscopic super- exchange 
paths between the chains (running along the c-direction) 
in CuGeOa. The three-dimensional situation is depicted in 
Ref. @. 

To determine the remaining parameters we equate four 
different experimental points with the corresponding pa- 
rameter dependent dispersion values given by Eq. (|). 
The parameters are fixed by solving the resulting sys- 
tem of equations. For (3 = 0.3 and —0.3 Figs. ||, || show 
the resulting dispersion curves in c*- and in b*-direction, 
respectively, using all hj <n calculated. 
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FIG. 6. Dispersion ui(k a = 0,kb,k c = 0) in b*- direct ion. 
Otherwise as in Fig. ^ 

As can be seen from Figs. ^| and ^| the plain series up 
to 10 th order provides excellent fits. Yet one realizes that 
the parameter values still change on passing from order 
to order. So it appears that even at 10 th order the results 
are not quantitative. In order to obtain quantitative re- 
liable results we adopt a systematic extrapolation in the 
order. In each order I € {3,4, ...10} we determine the 
optimum fit parameters. For illustration, Fig. [7] shows 
results for (3 = 0.3. 
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values for \i are listed in Tab. III. 
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FIG. 7. Dependence of the parameter values on the per- 
turbation order I at j3 = 0.3. The quasi-constant behavior of 
fj, is found for all /3- values checked. The lines are fits to the 
data according to Eq. (^). The solid lines consider all points, 
the dashed ones only the last four points. 

Assuming exponential convergence we use 

/(/) = X - be~ 2cl , (4) 

where X is the asymptotic value of the parameter con- 
sidered and b and c are constants. The choice (|]) is 
motivated on one hand by its obvious applicability (see 
Fig. [?]). On the other it stems from the fact that 
CuGeC>3 is a quasi one-dimensional gapped spin sys- 
tem. So one expects the magnetic correlations to drop 
exponentially with distance. Furthermore the order I 
determines the maximum distance over which correla- 
tions occur (cf. Ref pd] , namely I counted in dimers 
or 21 counted in spin sites. Hence the constant c in 
Eq. (Q) can be understood as the inverse of a correla- 
tion length £. With the usual relation £ « t>s/A for 
one dimensional systems we obtain c ~ 1/6 based on 
the rough estimates vs = 7r/2 • J(l — 1.12ao) JlTj and 
ao ~ 0.3; J « 12meV; A w 2meV. This is indeed what is 
found (cf. Tab. I, II) so that we judge our extrapolations 
as being well justified. 

In Fig. |?] the extrapolations are depicted by lines. The 
solid lines were obtained by using all calculated param- 
eter values. The dashed lines are obtained from on the 
last four points, i.e. the results in order 7, 8, 9, and 10. 
The deviation between these two extrapolations are used 
as a measure for the extrapolation error. This procedure 
is carried out for a, A, and J. 

There is no systematic dependence of [i on the order 
I. The parameter fj, oscillates between the two thin hor- 
izontal lines in Fig. [?]. So we take the average of these 
two bounds as our estimate for /x and their difference as 
the error in the determination of fi. 

Tabs. I and II summarize the results of the fits for the 
parameters a, A and J for four preset values of (3. The 





All points considered 


last four points considered 




X \ b \ c 


X \ b \ c 


(3 = 0.3 


a 


0.245 


2.61 


0.249 


0.236 


8.24 


0.338 


A 


0.867 


0.501 


0.155 


0.840 


8.66 


0.418 


J/moV 


13.6 


11.0 


0.164 


13.2 


105 


0.357 


P = 0.22 


a 


0.232 


3.27 


0.268 


0.228 


9.50 


0.343 


A 


0.863 


0.681 


0.184 


0.842 


13.6 


0.444 


J/moV 


13.1 


11.3 


0.171 


12.8 


68.9 


0.324 


(3 = 


a 


0.218 


4.16 


0.294 


0.226 


9.20 


0.334 


A 


0.862 


0.902 


0.213 


0.848 


14.94 


0.442 


J/mcV 


12.8 


11.1 


0.174 


12.6 


44.48 


0.287 


(3 = -0.3 


a 


0.212 


4.48 


0.300 


0.222 


9.94 


0.337 


A 


0.863 


0.974 


0.218 


0.849 


16.5 


0.448 


J/moV 


12.7 


10.9 


0.173 


12.5 


42.2 


0.282 



TABLE I. Extrapolated parameter values X according to 
Eq. (^). The experimental points we used in the fit process 
for this table are (cf. Figs ^, ^) [(kb, fc c );aj(k)/meV]: 
[(0,0); 2.1], [(0,0.05); 4.55], [(0, 0.25); 15.7], [(1, 0); 5.78]. 





All points considered 


last four points considered 




X | 6 | c 


X | b | c 


P = 0.3 


a 


0.297 


1.41 


0.215 


0.309 


1.18 


0.187 


A 


0.868 


0.423 


0.197 


0.877 


0.404 


0.165 


J/meV 


14.3 


8.88 


0.142 


13.9 


37.4 


0.269 


P = 0.22 


a 


0.301 


1.41 


0.208 


0.318 


1.15 


0.175 


A 


0.867 


0.850 


0.274 


0.886 


0.526 


0.175 


J/meV 


13.6 


12.0 


0.191 


13.8 


15.2 


0.191 


p = o 


a 


0.308 


1.34 


0.191 


0.323 


1.30 


0.173 


A 


0.900 


0.38 


0.140 


0.896 


0.754 


0.191 


J/meV 


13.6 


8.00 


0.133 


13.7 


9.31 


0.139 


P = -0.3 


a 


0.314 


1.28 


0.180 


0.326 


1.37 


0.173 


A 


0.913 


0.369 


0.127 


0.903 


0.832 


0.192 


J/meV 


13.6 


7.69 


0.127 


13.6 


9.12 


0.136 



TABLE II. Same as in Tab. I based on different exper- 
imental points: [(0,0); 2.1], [(0, 0.05); 4.35], [(0, 0.25); 15.7], 
[(1,0); 5.78]. 

A closer inspection of Figs. |^ and || reveals that we 
are confronted with a certain arbitrariness of which fit 
we should favor. The experimental errors enhance this 
problem. The filled circles in the range of small wave vec- 
tors in Fig. ^| represent experimental points which have 
been measured with a high degree of precision. Thus it is 
reasonable to fit the theoretical curve as well as possible 
to these points. Fig. || shows an enlargement of this re- 
gion. The solid line is the 10 th order fit result for (3 = 0.3 
as the solid line in Fig. ^. The depicted arrow indicates 
the experimental point (k c = 0.05, uo = 4.55meV) used to 
obtain Tab. I. 
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exp. points 
exp. points 
exp. points 

fitting-point: (0.05, 4.55 meV) 
fitting-point: (0.05, 4.35 meV) 



0.10 



k c /[r.l.u.] 



FIG. 8. Enlargement of Fig. ^ for small wave vectors at 
(3 — 0.3. The solid curve is the same as in Fig. leading to 
Tab. I. The dashed line corresponds to a 10 th order fit where 
the point (k c = 0.05, u) = 4.55meV) indicated by the arrow 
is replaced by (k c = 0.05, u — 4.35meV), leading to Tab. II. 
The filled circles correspond to highly accurate experimental 
points for which the error bars are of the size of the symbols. 

A likewise well suited curve, however, is produced if 
one uses the point (k c = 0.05, w = 4.35meV) for the fit 
keeping the other points (Tab. II). It is not possible to 
prefer one of the two lines in Fig. || to the other on the 
basis of their agreement to the experimental data. Hence 
we choose these two fits as the bounds within which all 
fits are acceptable. The corresponding parameter values 
Xi (fit 1) and X 2 (fit 2) provide an interval [Xl,A 2 ] 
which we expect to contain the true model parameter X. 
Hence the latter is estimated by 



X = (X 1 + X 2 )/2 ± AX 



with A, = 1/2 (X, 



all points 



■X 



last 4 pointsx 



AX,; = \Xi 



all points i 



and AX = max{|X - X x \, AX U AX 2 }. The 



X, 

results are summarized in Tab. III. 



parameter 


interval 


parameter 


interval 


= 0.3 


a 


0.27(4) 


a 


0.25(3) 


A 


0.86(2) 


S 


0.08(1) 


M 


0.27(1) 


Mo 


0.29(1) 


J/meV 


13.8(5) 


Jo /meV 


12.8(6) 


(3 = 0.22 


a 


0.27(4) 


a 


0.25(3) 


A 


0.86(2) 


S 


0.08(1) 


M 


0.21(1) 




0.23(1) 


J 1 me V 


13.4(4) 


Jo /meV 


12.5(5) 


(3 = 0.0 


a 


0.27(5) 


a 


0.25(5) 


A 


0.88(3) 


S 


0.07(2) 


A* 


0.13(1) 


Mo 


0.14(1) 


J/meV 


13.2(6) 


Jo /meV 


12.4(7) 


(3 = -0.3 


a 


0.27(5) 


a 


0.25(5) 


A 


0.88(3) 


S 


0.06(2) 


A 1 


0.08(2) 


Mo 


0.09(1) 


J/meV 


13.1(5) 


Jo /meV 


12.3(7) 



TABLE III. Final parameter intervals resulting from Tabs. I 
and II for three different values of /3. 

For the readers' convenience Tab. Ill also gives the 
results in the more commonly used parameters S, cvo, Mo> 
Jo and p. This notation is connected to the one used so 
far in this article by 



J 



Jo(l + <5), 

O'O 



a 



1 



A: 

M 



1 



1 + 5 
Mo 



1 



(5) 



It corresponds to the Hamiltonian depicted in Fig. ^. 
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FIG. 9. Alternative notation for the couplings in the 
dimerized phase of CuGeOs. The couplings are denoted rela- 
tive to the average nearest-neighbor coupling Jo in the chains. 
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IV. IMPLICATIONS FOR THE SUSCEPTIBILITY 

The temperature dependence of the homogeneous sus- 
ceptibility x(T) is often used to determine the parameters 
of CuGe0 3 . Already the Curie- Weiss tempera- 

ture provides valuable information on the sum of the 
coupling constants. This is particularly useful to detect 
frustration. The dispersions are governed by the differ- 
ence of the direct and the frustrating coupling whereas 
x(T) at larger temperatures is more influenced by the 
sum of direct and frustrating coupling. 

The analysis of the Curie- Weiss temperature alone 
bears some risks. It is easy to calculate but difficult to de- 
termine experimentally since it has to be deduced from 
values at high temperatures where x(T) is fairly small 
and hence strongly influenced by background effects (van 
Vleck, diamagnetism) or by slight structural changes. 

A convincing fit for temperatures above 50K is given 
by Fabricius et al. in Ref. |2(| on the basis of frus- 
trated chains. The inclusion of interchain couplings, 
however, would spoil the excellent agreement and a re- 
determination of the constant would be necessary. A 
description of x{T) on the basis of a two-dimensional 
model has not been done except for a consideration of 
the two leading coefficients in an expansion in 1/T in [I]. 
In Fig. [l^ we show the same high quality experimental 
data as in Ref. Q and compare it to theoretical curves at 
four values of j3. The theoretical curves are obtained by 
computing a [4,5] Dlog Pade approximant xq(T) based on 
the high temperature series provided in Ref. p2| for the 
frustrated chains. This procedure provides excellent re- 
sults down to T « J/5 [E3[ . The asymptotic behavior of 
the approximant is chosen such that Xo(T) vanishes lin- 
early on T — > as is to be expected for a two-dimensional 
massless antiferromagnet. Besides this feature the two- 
dimensionality is incorporated on a chain-mean-field level 



0.0020 



X(T) = 



Xo(T) 



l + 2 Mo (l + 2/3)xo(T) 



(6) 



This relation is exact in linear order in /i. Estimates of 
corrections quadratic in /i indicate that they are negligi- 
ble for the values of \x and p in which we are interested. 




0.0000 



400 600 

T/[K] 



1000 



FIG. 10. Comparison of experimental data in b-direction 
po[ and theoretical susceptibilities for various values of the 
interchain frustration. The g value used is 2.26 |2~i|| 

From the results in Fig. |l^ it is evident that the inter- 
chain frustration cannot be neglected. Only for a finite 
value of about 0.22 a very good agreement can be ob- 
tained. The agreement to the frustrated chain model 
p(if is still better since the position of the maximum is 
also reproduced. But on the basis of the neutron scat- 
tering results (lj] it is undoubtful that CuGeOa is a two- 
dimensional substance. Furthermore, it must be consid- 
ered that the previous fit [po[ was a two-parameter (J 
and ao) fit whereas only one parameter is fitted to 
obtain Fig. [l(]. The other parameters (Jo,aojMo) were 
determined from an entirely different experiment. Hence 
the agreement for x(T) corroborates also the validity of 
the parameters determined in the preceding section. 



V. DISCUSSION 

We will first discuss our results and propose a set of pa- 
rameters. Then we will put these results into the context 
of other results in the literature. 

Let us consider the remaining difference between ex- 
periment and theory concerning x{T) in Fig. [l^. There 
are four conceivable sources for it. The first are experi- 
mental inaccurcies. We are not in the position to judge 
this aspect. We just like to comment that from the re- 
sults and error bars in Fig. || it is obvious that the ex- 
perimental data is not completely consistent so that this 
explanation is possible. 

Second, it is conceivable that the couplings change 
across the transition, i.e. the intrachain frustration at 
T m (where the dispersions are measured) is differ- 
ent from the one above T SP (where x{T) is measured). 
Since we are considering a spin-Peierls transition there is 
definitely a structural change. So far, however, the as- 
sumption that only the dimerization changes worked well. 
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The structural changes in the transition are very small 
iful whereas the changes needed to explain the discrep- 
ancy are of the order of 20 to 30 % (assuming a change 
in the intrachain frustration). Estimates point into the 
direction that the changes on the couplings are unimpor- 
tant pq| . Yet the estimates concern in the first place 
the nearest-neighbor coupling only. Quantitative ab ini- 
tio calculations of the frustration are very difficult jlB) , 
even more so for changes of the frustration. So, again, 
this explanation is perhaps not the most promising but 
cannot be excluded either. 

Third, the influence of the phonon dynamics is to be 
considered, ft is shown that spin-Peierls systems can be 
unitarily transformed in such a way that an effective spin 
model remains at low energies uncoupled to the phonons 
p?j] . The effective couplings in a single chain model be- 
come temperature dependent so that this may account 
for the discrepancy. But it turns out for non-resonant 
phonons (ui > J) that these effects leave the susceptibil- 
ity fairly unchanged. This is so since these effects become 
significant for relatively large temperatures where x{T) 
depends only on the sum of all couplings which is un- 
changed by the unitary transformation (this is observed 
in (VO)2P207 ^8|). So this reason appears rather un- 
likely even though it looked plausible at first sight. 

Fourth, one has to think about any kind of precursors 
of the spin-Peierls transition. By this we mean on one 
side the critical fluctuations which appear in a narrow 
region (ss 3K) around the spin-Peierls temperature p6| . 
On the other side, we mean any precursor which goes 
beyond a purely static spin model. Experimentally, a 
finite lattice correlation length can be detected already 
at T » 40K far above the actual transition p9| . From 
there on deviations from the behavior of a static spin 
model should be observable. In the adiabatic limit, for 
instance, the fluctuation yield already a reduction of the 
susceptibility J3(]]. What happens in the antiadiabatic 
limit has not yet been investigated for a two or higher 
dimensional model. The mapping in Ref. J27J leads to 
four-point interchain couplings the influence of which is 
unclear so far. 

In view of the above mentioned possible pitfalls of the 
static model the agreement in Fig. |fC| is already very con- 
vincing. Summarizing our results we propose the param- 
eters given for (3 = 0.22 in Tab. Ill to be the ones deduced 
from the dispersion data. Assessing the reliability of our 
estimates, we redo the analysis of the susceptibility for 
cto = 0.28 (the upper bound of our estimate for ag) with 
the corresponding value of Jo = f 2.8meV. Then the op- 
timum x(T) is obtained for j3 — 0.15; the x(T) curve is 
almost identical to the one shown in Fig. [n]. So the value 
of [3 cannot be determined very precisely, but it should 
be in the range j3 = 0.2(1). A certain dicrepancy between 
the optimum parameters for the T = dispersion data 
and for the x{T > Tsp) data remains. 

We split the comparison of our findings to previous 
works into three groups. The first comprises the analy- 
ses on the basis of a one- dimensional model 



The most striking difference to the results for static spin 
models |T^- 12(]] is that the dimerization 6 is not of the 
order of 1% but significantly larger. This is not astound- 
ing since it has been noted already in [I] that the gap 
is lowered by the interchain hopping. Hence the neglect 
of the latter requires to lower the gap otherwise, i.e. by 
assuming a lower dimerization. 

Our intrachain frustration is slightly larger than the 
one of Castilla et al. (0.24), but significantly smaller than 
the value of Riera and Dobry (0.36) or the value of Fabri- 
cius et al. (0.35). Fabricius et al. showed that the value 
cto = 0.24 is too small for a single chain model. The dif- 
ference between the larger frustration value in the single 
chain model to our value results directly from the inter- 
chain coupling. As can be nicely seen in Eq. (^|) the inter- 
chain coupling lowers the susceptibility without changing 
(in the chain mean-field approximation) the position of 
its maximum. The one dimensional models favor a larger 
intrachain frustration and a concomitant larger coupling 
J in order to reduce the magnitude of the susceptibility. 

The claim by Wellein and coworkers that the dimer- 
ization experimentally found to be larger than would fit 
to a static ID model |32| is due to the phonon dynamics 
is not compelling. They use a root-mean-square defini- 
tion of the dimerization which naturally provides larger 
values for the dimerization since it includes all the fluc- 
tuations. The dispersion perpendicular to the chains, 
however, is an unambiguous experimental fact. Further- 
more, Trebst and coworkers [ |33"| do not find a substantial 
gap renormalization for parameters relevant for CuGe03 
even though one should take care of different schemes to 
couple the phonons. 

Let us turn to the second group of papers consider- 
ing the essentially two-dimensional character of CuGcOs. 
The first work used a bond-operator technique fl34| . No 
frustration was considered, hence rather small values of 
J = I0.2meV and a rather small interchain hopping 
(i Ri 0.06 resulted. The same technique was also ap- 
plied later again by Brenig p5|] including frustration. It 
turned out, however, that only A(l — 2a) and /i(l — 2[3) 
matter on the free-boson level. Hence an independent 
determination of the frustration is not possible. Using 
additional input (S = 0.012) the values cto — 0.059 and 
/i(l — 2/3) = 0.054 were obtained. In view of the exten- 
sive comparisons to numerical results made in Ref. |35| ] 
it appears that the bond-operator method overestimates 
the influence of additional couplings such as dimerization 
or frustration. Generally, the values for dimerization or 
frustration tend to be too low. This is confirmed by our 
findings in the present work. 

Compared to [I] (ao = 0, J = 9.8mcV, S = 0.12, ^ = 
0.34, P = 0.3) the extended series on which our present 
analysis is based gives a much better handle on intrachain 
frustration, see Fig. Only in the present work, we are 
able to assess its value reliably. With respect to the inter- 
chain frustration, the present results agree qualitatively 
with those in [I] where such a frustration was proposed 
for CuGeOs first. The use of susceptibility information 
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has been improved in the present work since the whole 
x(T) curve is used, not only the leading coefficients. 

Bouzerar et al. have carried out an estimate leading 
to results not too far from ours: 6 = 0.065, ao = 0.2, J = 
12.2mcV, /i = 0.15. They used just linear order in the in- 
terchain hopping without interchain frustration and some 
square-root averaging with numerical results for chains 
to describe the dispersion. The intrachain frustration 
(ao = 0.2) could only be taken from the Curie-Weiss 
constant. The resulting x(T) nas similarities with the 
experimental one. 

The third group comprises ab-initio calculations of the 
exchange couplings and of the spin-phonon couplings. 
Microscopic calculations |36|,|l6| find relatively large val- 
ues of the dimerization between 0.07 and 0.2 in agreement 
with our findings. (Even though there is also a different 
result [j37| ) Very important for our work are recent results 
by Drechsler and coworkers supporting the existence of 
interchain frustration jlB). Werner and coworkers f2(f | 
estimate a large dimerization from the spin-phonon cou- 
plings and the shifts of the ions (<5 = 0.11). From the bal- 
ance of elastic and magnetic energy in the D phase they 
obtain without frustration a lower bound of S > 0.044. 
Assuming critical frustration ao = 0.2412 they find even 
S > 0.078 which fits very nicely to our findings. 

In summary, we provide by the present work a de- 
termination in great detail of the coupling parameters 
(f3 = 0.2(1) and right column in Tab. Ill under (3 = 0.22) 
of CuGeOs based on a static dimerized spin model at 
T = 0. The experimental input comes from inelastic 
neutron scattering. The implications of the parameters 
found for the susceptibility are also studied. Very good 
agreement could be obtained fitting the interchain frus- 
tration appropriately. A small discrepancy at low tem- 
peratures around 50K indicates that the static spin model 
is probably insufficient to describe CuGeC-3 completely. 
By this work, we proved the outstanding possibilities of 
high-order series expansions (around the dimer limit or 
around the limit T = oo) in the analysis of experimental 
data. 
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